Signatures of Inelastic Scattering in Coulomb-Blockade Quantum Dots 
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We calculate the finite-temperature conductance peak-height distributions in Coulomb-blockade 
quantum dots in the limit where the inelastic scattering rate in the dot is large compared with the 
mean elastic tunneling rate. The relative reduction of the standard deviation of the peak-height 
distribution by a time-reversal symmetry-breaking magnetic field, which is essentially temperature- 
independent in the elastic limit, is enhanced by the inclusion of inelastic scattering at finite tem- 
perature. We suggest this quantity as an independent experimental probe for inelastic scattering in 
closed dots. 



Quantum coherence of the electrons in ballistic quan- 
tum dots leads to distinct signatures in the mesoscopic 
fluctuations of their conductance 0). These signatures 
are affected by decoherence and inelastic scattering of 
the electrons at finite temperature. 

Various methods can be used to measure the dephas- 
ing rate of the electrons in open dots H| . In particular, 
the weak- localization effect (i.e., the reduction of the av- 
erage conductance in the absence of a magnetic field) is 
sensitive to decoherence, and the suppression of this ef- 
fect has been used to extract from experimental data the 
temperature dependence of the dephasing time. 

In almost-isolated dots, the dephasing and inelastic 
scattering rates are expected to vanish at low excitations 
0, 0j. It is, however, more difficult to determine experi- 
mentally the dephasing rate in such dots. An effect anal- 
ogous to weak localization was predicted in the average 
Coulomb-blockade peak height of a closed (i.e., almost- 
isolated) chaotic dot fij, and, in the limit of pure elas- 
tic scattering, was found to be essentially independent 
of temperature 0. A recent experiment measured the 
weak-localization effect in closed dots and found good 
agreement with theory for temperatures kT below the 
mean level spacing A ||. However, this effect was found 
to be suppressed for kT > A. Recently, the opposite limit 
where the mean inelastic scattering rate T m is large com- 
pared with the mean elastic rate r c i [| was studied, and 
the weak-localization effect in the average peak heights 
was found to be suppressed with increasing temperature 
p0[ . This suppression was stronger than what was ob- 
served in the experiment, and it was concluded that the 
experimental data are consistent with Ti n < r c i <C A. 
Ref. has generalized the master-equations approach 
of Ref. M to study the effect of inelastic scattering on the 
average peak height in the crossover regime Ti n ~ r e i- 
This allows one to determine from the data a mean in- 
clastic scattering rate. 

Here we study additional signatures of inelastic scat- 
tering on the conductance peak-height statistics in closed 
Coulomb-blockade quantum dots. We calculate the peak- 
height distributions at finite temperature in the inelastic 
limit r in » r c i, and compare them with the correspond- 
ing distributions in the elastic limit. In particular, we 



find that the decrease in the standard deviation of the 
conductance peak height upon application of a magnetic 
field (measured in units of the peak-height standard de- 
viation in the presence of a magnetic field) is sensitive to 
inelastic scattering. This reduction in the standard devi- 
ation is temperature-independent in the elastic limit but 
is enhanced as a function of temperature in the inelastic 
limit, and therefore can be used as an experimental sig- 
nature of inelastic scattering that is independent of the 
weak-localization suppression. 
In Ref. U, 

the ratio of the standard deviation of the 
peak height to its mean value was measured as a function 
of temperature in the presence of a magnetic field and 
compared with theory. It was found that the experimen- 
tal values are suppressed in comparison with theoretical 
values calculated in the elastic limit [[l3) , and the suppres- 
sion is stronger at higher temperatures. To determine 
whether this suppression can be explained by inelastic 
scattering, we calculate the ratio of standard deviation 
to mean peak height in the inelastic limit. While we find 
some suppression (compared with the elastic limit), it is 
relatively small and insufficient to explain the observed 
data. 

We consider a closed dot that is weakly coupled to the 
leads, where the average width r o i for elastic decay of 
an electron into the leads is much smaller than the mean 
spacing A of the single-particle levels in the dot. We also 
assume that the thermal energy kT is much smaller than 
the charging energy e 2 /C, where C is the capacitance 
of the dot (Coulomb-blockade regime). To identify the 
effect of inelastic scattering of the electrons within the 
dot on the linear conductance, we compare two limits. 
In the elastic limit, the inelastic decay rate is assumed to 
be negligible compared with the elastic rate Ti n <C r e i, 
and an electron decays into one of the leads prior to any 
inelastic scattering event inside the dot. In the opposite 
limit of strong inelastic scattering, T m 3> r e i, an electron 
tunneling into the dot has sufficient time to thermalize 
because of inelastic processes before leaving the dot. 

Model. We further assume the limit r c i,ri n < kT of 
sequential tunneling that can be described by a master- 
equations approach In the Coulomb-blockade regime, 
the conductance displays sharp peaks as a function of 
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gate voltage. The conductance G around a Coulomb- 
blockade peak can be written as 



g(T,E F ) 



(1) 



where g{T, Ep) is a dimensionless conductance at temper- 
ature kT and effective Fermi energy E-p (the latter can be 
tuned by varying the gate voltage). For fcT, A <§; e 2 /C 
and in the elastic limit, g — g c \ can be expressed as a 
sum of contributions from several levels A in the dot || 
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where T A is the partial elastic width of the level A to 
decay to the left (right) lead, and T e i = rl 



T' A . The 



contributions of the levels A in Eq. 
factors 



are weighted by 



w x = 4f(AF N - E F )(n x ) N [l - f{E x - Ep)] , (3) 

which are a function of both T and E-p JT3| . Here, / is the 
Fermi-Dirac distribution function, AFn = Fn — -Fjv-i, 
where Fn is the canonical free energy of N noninter- 
acting particles, and (n\)N are the canonical occupation 
numbers of the levels A in the dot containing N electrons. 

On the other hand, in the inelastic limit Fj n 3> r o i, 
there is full thermalization of the electrons among the 
levels in the dot and g ~ g m with H 



2 (EA^r A )(E^r; 
r cl E,^(ri + r r j 



(4) 



Note that the inelastic conductance (|j) does not depend 
explicitly on r m (in the inelastic limit). 

In a chaotic dot (e.g., with irregular shape), the single- 
particle levels and wave functions follow random-matrix 
theory (RMT) statistics (I). In particular, the partial 
elastic decay widths are fluctuating quantities that are in- 
dependent and distributed according to a Porter-Thomas 
law P PT (r) oc r /3 / 2 - 1 exp(-/3r/r el ), where the symme- 
try parameter /? indicates the presence (/3 = 1) or ab- 
sence (/3 = 2) of time-reversal symmetry. We assume 
that all levels have a common average partial width, i.e., 

iFrT \(r) 

Iv = Tj2 for any pair of levels A and /i, and that the 
leads are symmetric: T\ = T l x . 

Method. We calculate the conductance peak-height 
distributions in the inelastic limit and compare them with 
the corresponding distributions in the elastic limit. In an 
RMT approach, we model the single-particle spectrum 
{E\} of the quantum dot by the spectrum of a matrix 
that belongs to the corresponding Gaussian ensemble, 
i.e., the Gaussian orthogonal ensemble (GOE) for con- 
served time-reversal symmetry (/? = 1), and the Gaus- 
sian unitary ensemble (GUE) for broken time-reversal 
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FIG. 1: Left panels: distribution -P(<?max) of the dimensionless 
peak height g max (see Eq. (Q) in the elastic (dashed lines) 
and inelastic (solid lines) limits at temperatures kT = 0.5A 
(top panel) and kT — A (bottom panel) . Right panels: mean 
g ma x (top) and standard deviation <r((7 max ) (bottom) of the 
peak height j max in the elastic (dashed lines) and inelastic 
(solid lines) limits. 



symmetry {(3 — 2). The widths T r x } are taken as the 
squares of corresponding eigenvector components of the 
random matrix. For a given number N of electrons, both 
the canonical free energy Fn and the occupation num- 
bers (n\) n are calculated exactly using particle- number 
projection Jig]] (projecting from grand canonical parti- 
tion functions with imaginary chemical potential). The 
elastic and inelastic conductance peaks at temperature 
T are then calculated as a function of the effective Fermi 
energy E-p from Eqs. (J2j) and (Q), respectively. Finally, 
the height <? max of the corresponding conductance peak at 
temperature T is calculated by maximizing the function 
g(T,Ep) with respect to Ep. The (numerical) distribu- 
tion P(<7max) of the conductance peak heights is obtained 
by calculating the peak heights for an ensemble of RMT 
spectra and corresponding widths. 

Results. We have carried out calculations for both non- 
degenerate and spin-degenerate single-particle spectra. 
In both cases A is the mean-level spacing per single- 
electron state (and is therefore half the orbital spacing 
in the spin-degenerate case) . In the following, we assume 
non-degenerate energy levels unless otherwise stated. 

Fig. |l] illustrates the temperature dependence of the 
distributions -P(gmax) of conductance peak heights g m &x 
in both the elastic and inelastic limits. In the left pan- 
els, the distributions -P(<? m ax) are shown at temperatures 
of kT = 0.5A and A. In the right panels, the mean 
conductance peak height g max and standard deviation 

" 2 ) are shown versus tem- 
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CT(5 m ax)(fT (ffmax) = 5: 

perature kT. The elastic quantities are shown by dashed 
lines and the inelastic quantities by solid lines. All the 
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results in Fig. |] are for the case of broken time-reversal 
symmetry (/? = 2), but similar behavior is found for the 
case of conserved time- reversal symmetry (/3 = 1). 

In the limit of kT <C A, only a single level A con- 
tributes to the various sums in Eqs. (||) and (||) and 
,9ci = ffin- In this limit, the distributions of the elastic 
and inelastic conductance coincide, and the mean peak 
height and its standard deviation take the analytically 
known GUE values for one- level conduction, g max = 1/3 
and o-(sw) = 2^5/15 « 0.2981 (the GOE values are 
1/4 and cr(g max ) = V2/4 « 0.3536). At finite 



and the standard dc- 



<?max 

temperatures, both the mean ~cf^ 
viation a(g ma , x ) are found to be enhanced by inclusion of 
inelastic scattering. 
In recent work 



10 1, the finite-temperature suppres- 



sion of the weak-localization effect in the average con- 
ductance peak height was suggested as a signature of 
inelastic scattering in the dot. Quantitatively, the weak- 
localization effect can be described by the increase of the 
average peak height upon breaking of time-reversal sym- 
metry (relative to the GUE average value) 
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^max 
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-GUE 



(5) 



At low temperatures kT < A, a = 0.25 || In the 
elastic limit, a is fairly independent of temperature |7j, 
as is shown by the open symbols in the lower panel of 
Fig. H The circles correspond to a non-degenerate single- 
particle spectrum and the diamonds to a spin-degenerate 
spectrum. As the temperature increases, more levels 
contribute to the elastic conductance but the ratio be- 
tween the GOE and GUE averages of the peak height 
remains approximately constant, since the average over 
several contributing levels commutes with the ensemble 
average. A small suppression in the temperature range 
kT/A < 0.8 is due to the suppressed occurrence of close 
lying levels in the GUE CclSC, ctS pointed out already in 
Ref . jllj] . Because level repulsion is stronger in the GUE 
than the GOE, the contribution of an excited level in the 
dot to the elastic conductance is smaller on average for 
the GUE, and a decreases somewhat. 

On the other hand, inelastic scattering is found to re- 
duce a at finite temperature and a — > for kT 3> A. 
As a reminder, we show in the lower panel of Fig [2] the 
temperature dependence of a in the inelastic limit (solid 
symbols). While the data shown in Ref. |HJ was ob- 
tained assuming equidistant levels, our results are calcu- 
lated from RMT spectra yielding a stronger suppression 
of the inelastic a at low temperatures. 

Here we suggest an independent signature of inelastic 
scattering that is experimentally accessible. We define 7 
as the reduction of the standard deviation of the conduc- 
tance peak height (relative to the GUE standard devia- 
tion) when a time-reversal symmetry-breaking magnetic 




FIG. 2: The quantities 7 (upper panel, Eq. (H)) and a (lower 
panel, Eq. (B)) versus temperature kT in the elastic (open cir- 
cles) and inelastic (solid circles) limits. The diamonds are the 
corresponding results for a spin-degenerate spectrum. Notice 
the sensitivity of both 7 and a to inelastic scattering. 



field is applied, 
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In the upper panel of Fig. we show 7 versus temper- 
ature for both the elastic (open symbols) and inelastic 
(solid symbols) limits. For kT < A, 7 = 3\/5/(4V2) - 
1 rs 0.1859, characteristic for conduction involving a sin- 
gle level in the dot. In the elastic limit, 7 is essentially 
independent of temperature. However, at finite tem- 
perature, 7 is sensitive to inelastic scattering. In the 
inelastic limit, 7 increases from 0.1859 at kT <C A to 
ss 0.42 ± 0.01 at kT > A. The high temperature limit 
is obtained using the approximation w\ « uk /2, where 

w a°' = V cosh 2 [(i? a — Ep)/2kT] are the thermal weights 
in the absence of charging energy 

The quantity 7 can be calculated in the crossover 
regime r in ~ r c i using the generalized master-equations 
approach of Ref. O]. It can therefore be used as an 
independent experimental probe from which the finite- 
temperature inelastic scattering rate can be determined. 

Another dimensionless quantity of relevance to the ex- 
periments is the ratio between the standard deviation 
c(gmax) and average <7max of the peak heights. This ratio 
was measured versus temperature in Ref. |l^] and found 
to be significantly suppressed at higher temperatures in 
comparison with RMT calculations jl3| . However, these 
predictions were based on the assumption of pure elas- 
tic scattering, and it was conjectured that the experi- 
mentally observed suppression might be due to inelastic 
scattering. Fig. |^ shows the temperature dependence of 
"■(SmaxJ/ffmax for both the elastic (dashed lines) and in- 
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FIG. 3: The ratio o"(g m ax) /<7max between the standard devia- 
tion and average value of the peak height versus temperature 
kT in the elastic (dashed lines) and inelastic (solid lines) limits 
for both GOE (left) and GUE (right) symmetries. The sym- 
bols in the right panel are the experimental data of Ref. 1l2j . 



elastic (solid lines) limits and for both conserved (left 
panel) and broken (right panel) time-reversal symmetry. 
The symbols in the right panel are the experimental re- 
sults jp|. While inelastic scattering enhances the sup- 
pression of the ratio cr(<7max)/<7max at finite temperature, 
it still cannot explain the experimental findings at the 
present stage. We note, however, that even at low tem- 
peratures there is a discrepancy between the experimen- 
tal data and the analytically known value for one-level 
conduction {a(g max )/g^ = 2V5/5 » 0.8944 for (3 = 2). 
As a probe of inelastic scattering, the parameters a and 
7 turn out to be far more sensitive than the ratio between 
the standard deviation and mean value. The increases of 



scattering inside the dot on the distribution of the con- 
ductance peak heights in the Coulomb-blockade regime. 
We found that the quantity 7 measuring the decrease in 
the standard deviation of the conductance peak height 
upon breaking of time-reversal symmetry (relative to the 
value of the standard deviation for broken time-reversal 
symmetry) is sensitive to inelastic scattering of electrons 
within the dot. This quantity can serve as an indepen- 
dent experimental probe of inelastic scattering, in addi- 
tion to the quantity a that was suggested and measured 
recently. 
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panels of Fig. cancel out to a large extent in the quan- 
tity <7(.9max)/5max, resulting in a reduced sensitivity of 
the latter to the inelastic scattering rate. 

Conclusion. We have investigated the effect of inelastic 
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